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ABSTRACT 

We report the Swift discovery of the nearby long, soft gamma-ray burst GRB 100316D, 
and the subsequent unveihng of its low redshift host galaxy and associated supernova. 
We derive the redshift of the event to be 2; = 0.0591 ± 0.0001 and provide accurate 
astrometry for the GRB-SN. We study the extremely unusual prompt emission with 
time-resolved 7-ray to X-ray spectroscopy, and find that the spectrum is best modelled 
with a thermal component in addition to a synchrotron emission component with a 
low peak energy. The X-ray light curve has a remarkably shallow decay out to at 
least 800 s. The host is a bright, blue galaxy with a highly disturbed morphology and 
we use Gemini South, VLT and HST observations to measure some of the basic host 
galaxy properties. We compare and contrast the X-ray emission and host galaxy of 
GRB 100316D to a subsample of GRB-SNe. GRB 100316D is unlike the majority of 
GRB-SNe in its X-ray evolution, but resembles rather GRB 060218, and we find that 
these two events have remarkably similar high energy prompt emission properties. 
Comparison of the host galaxies of GRB-SNe demonstrates, however, that there is a 
great diversity in the environments in which GRB-SNe can be found. GRB 100316D is 
an important addition to the currently sparse sample of spectroscopically confirmed 
GRB-SNe, from which a better understanding of long GRB progenitors and the GRB- 
SN connection can be gleaned. 

Key words: gamma-ray burst: individual: GRB 100316D; supernovae: individual: 
SN2010bh 



1 INTRODUCTION 

The connection between long-duration gamma-ray bursts 
(GRBs) and Type Ic core-collapse supernovae (SNe) 
has long been established, beginning wi th the asso- 
ciation of GRB 980425 with SN 1998bw l|Galama et all 
1 19981 : iPian et al.l [l998). Subsequent associations between 
nearby GRBs and spectr oscopically-confirme d SNe include 
GRB 030329/SN 2003dh (iHiorth et al.l l2003l: IStanek et all 
I2OO3I ). GRB 031203/SN 20031w (^ Malesani et al.l |2004|) "aiid 
most recently GRB 060218/SN 2006aj l|Campana et all 
l2006l : IPian et al.l bood). These are spectroscopically con- 
firmed examples of nearby GRB-SN associations (out to 
z = 0.17). The characteristic 'bump' of a supernova has been 
noted in the p hotometric data for a dozen or more GRBs out 
to z ^ 1 (e.g. 'Ze h. Klose fc Hartmann"2005': 'Ferrero et al.' 
20061: IPc Ua Vall e et all I2OO6: Wooslcy & Bloom 2006; 



Tanvir et al.ll20ld ). while the majority of GRBs lie at higher 



redshifts where such signatures would be impossible to de - 
tect ((z) = 2.2 IJakobsson et ahlliooel : iFvnbo et al.ll2009l ). 
The GRBs with SNe are therefore rare, but provide a cru- 
cial insight into the GRB-SN connection and the progenitors 
of long GRBs. 

The GRBs with spectroscopically confirmed SNe are 
generally underluminous and subenergetic in comparison to 
a typical long GRB, though GRB 030329 is a notable excep- 
tion to this. The prompt emission of these bur sts has a lower 
energy spectral peak than a typical GRB (|Kaneko et al.l 
l2007l 'l and they are suggested to have less relativistic out- 
flows, or be viewed more off-axis. It has been suggested that 
the observed nearby {z < 0.1), underluminous/subenergetic 
GRBs, including GRB 100316D, may be drawn from a dif- 
ferent population to the cosmological GRB sample, moti- 
vated by estimates of the local GRB rate se veral times 
greater than the rate of cosmological GRBs (e.g. I Cobb et"al] 
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I2OO6I : IChapman etal] l2007l : iLiang et al.1 12007| ). Radio ob- 
servations of GRBs suggest that mildly relativistic ejecta 
and often high luminosities are the properties that distin- 
guish the supernovae associated with GRBs from the non- 
relativistic ordinary Type I b/c supernovae (|Kulkarni et al.l 
ll998l : ISoderberg et al.|[200^ . 

GRB 060218 was a landmark long-duration, low- 
luminosity, soft spectral peak eve nt, detected by the 
Swtft satellite (iGehrels et al.l |2004| ) with unprecedented 
multiwavelength coverage of the prompt emission. The 
prompt spectrum was found to comprise both the non- 
thermal synchrotron emission ascribed to most GRBs 
and thought to originate in the collision of fast-movin g 
shells within the GRB jet (e.g. iRees fc MeszarosI liggil ). 
and a thermal component. The presence of this thermal 
component, together with its evolution, led to the sug- 
gestion that we were observing t he shock break o ut of 
the supernova for the first time (jCampana et al.l l2006l : 
IWaxman. Meszaros fc Gampanal |2007|). However, the non- 
thermal emission did not differ greatly from that of the 
X-ray flash class of soft GRBs and an ou tflow speed 
close to the speed of light could be inferred l|Toma et al.l 
[2007), alternatively suggesting this to be an extension 
of the typical GRB population and not requiring sig- 
nific antly slower ejecta or any special (off - a xis) geome- 
try llGhiseUini. Ghirlanda fc Tavecchid 120071 ). iToma et all 
((2003) speculate that the long-duration, low-luminosity 
events such as this one may point to a different central en- 
gine for these events compared with more ty pical GRBs: a 
neutron star rather th an a black hole (see also lMazzali et al] 
I2OO6I : iFan et aLllioiol ). 

Also of relevance to this discussion are the supernovae 
for which no GRB was detected, yet for which mildly 
relativistic ejecta pro vide a good explan ation for the ra- 
dio data: SN2007gr (|Paragi et al.1 I2OI0I ) and SN 2009bb 
l|Soderberg et al.l |201G| ). As GRB jets are thought to be 
highly coUimated, it is expected that many will go unde- 
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tecte d if the emission is directed away from our line of sight 
(e.g. iMazzah et all 120051 ) . and a fraction will also lie below 
current detector sensitivity limits. ISoderberg et^ (j2010l ) 
estimate the fraction of Type Ib/c SNe with central en- 
gines at about one per cent, consistent with the inferred 
rate of nearby GRBs. In contrast, there are nearby GRBs 
for which an accompanying supernova could be expected 
but none has been det ected to extremely deep limits. GRBs 
060505 and 060614 (iFvnbo et al.l l2006l : iDella Valle et al.1 
I2OO6I : iGal-Yam et al.l I2OO6I ) are the two best examples of 
this, highlighting the need for a greater understanding of 
the relationship between GRBs and SNe. 

We presen t the discovery of a new GRB-SN, 
GRB 100316 E) jStamatikos et ai.l I2OI0I') associated with 
SN2010bh JWiersema et al.1 I2OI0I : iBufano et al.l I2OI0I : 



IChornock et al.ll2010l ). This is an unusually long-duration, 
soft-spectrum GRB positioned on a nearby host galaxy. 
Many of the properties of this GRB appear unlike those 
of typical GRBs, resembling rather GRB 060218. In this pa- 
per we report accurate astrometry for the GRB-SN and the 
redshift of the underlying host galaxy. We examine the GRB 
prompt and afterglow emission observed with all 3 instru- 
ments on-board Swift, and the broad characteristics of the 
host galaxy as observed with the Gemini South telescope, 
the Very Large Telescope (VLT) and the Hubble Space Tele- 
scope (HST). We compare these properties with the GRB 
sample as a whole and with a subset of GRB- SNe, to un- 
derstand the origins of the GRB emission and further our 
knowledge of the GRB-SN connection. 



2 OBSERVATIONS AND ANALYSES 

2.1 Swift 

On 2010 March 16 at 12:44:50 UT (hereafter Tg), the 
Swift Burst Alert Telescope (BAT. iBarthelmv et al.1 12OO5I ') 
triggered on and slew ed immediately to GRB 100316D 
IStamatikos et al.llioiol '). This was an image trigger, since 
no rapid, substantial rise in the count rate occurred, but 
a possible low-level peak about 100 s long was noted 
in Istamatikos et al.l (|2010l '). The X-ray Telescope (XRT, 
iBurrows et al.l l2005h on-board Swift began observing the 
field at 12:47:08.1 UT, 137.7 s after the BAT trigger, be- 
ginning with 10 s of Windowed Timing (WT) mode settling 
data and followed by 593 s of WT pointing data. After a 
time gap due to observing constraints, data-taking recom- 
menced at To +33 ks with XRT in Photon Counting (PC) 
mode. The best X-ray position for this source is the as- 
trometrically corrected XRT position of RA, Dec (J2000) 
= 07h 10m 30.63s, -56d 15' 1 9.7", with an uncer tainty of 
3.7" (radius, at 90% confidence fstarlin g et al."201(f). Swift's 
Ultra Violet-Optical Telescope (UVOT, Rorning ct al. 200i) 
began observing with the white filter 148 s after the BAT 
trigger and continued observations in u band, thereafter cy- 
cling repeatedly through all its lenticular filters. 

The Swift data were processed and analysed with the 
Swtft tools version 3.5, released as part of HEASOFT 6.8, 
and the latest calibration data at 2010 March were used. 

XRT WT spectra were grouped with a minimum of 
20 counts in each bin, allowing the use of statistics, 
while for the PC spectrum 1 count per bin was required 



and Cash statistics used. We fit the X-ray spectra over the 
energy range 0.3-10 keV and the 7-ray spectra over the 
range 15-50 keV (as this is a soft event very few counts 
are present a bove 50 keV) using the spectral fitting pack- 
age XSPEC (|Arnaud|[l996l ). We used the X-ray absorp- 
tion models phabs and zpha bs for which photoi onization 
cross-sections are ta ken from Verner et al. ll 19961) an d So- 
lar abundances from IWilms. Allen fc McCravl (|2000h . The 
Galactic absorptio n in the X-ray band is fixed at 7x10^" 
cm~^ throughout (|Kalberla et al.l [20051 ). We tested for any 
differences in the spectral fitting results from use of different 
grade selections or different energy ranges (excluding 0.3-0.5 
keV and/or including 50-150 keV) and confirm that these 
have no effect. 

Optical/UV photometry was performed with the UVOT 
tool UVOTSOURCE . At the position of the SN reported by 
iLevan et all (|2010l ). background subtracted count rates were 
extracted using a circular source aperture of 2.5". Source- 
free background was taken from a region situated outside 
of the host galaxy. The count rate was aperture corrected 
to 5" in order to be compatible with the UVOT calibra- 
tion (|Poole et aLlliboa ). With UVOT we detect the host 
galaxy, but no new point sources. The 2>a upper limit on 
any optical transient emission was taken to be the back- 
ground subtracted count rate plus three times the count 
rate error and was converted in to magnitude using the zero 
points provided in iPoole et al.l l|2008i) . At this time, galaxy 
subtraction cannot be performed since there are no UVOT 
pre- or post-explosion observations. The host galaxy count 
rate was determined using an elliptical source region, which 
enclosed the whole galaxy, and a source-free background re- 
gion. Photometry of the host was performed on a sum of all 
data from ~ 1 x 10^ s after the trigger onwards for each filter. 
Since the ellipse is similar in area to a 5" circular aperture 
(semi-major and -minor axes of 7.5" and 5.5" respectively), 
we have taken the standard UVOT co-incidence loss correc- 
tion to be the best estimate of the co-incidence loss, and the 
standard UVOT zero points have been used to convert the 
count rate into magnitudes. 



2.2 Gemini, VLT and HST 

Our first optical observations of GRB 100316D were ob- 
tained with Gemini-South, starting on March 16 23:53 UT. 
In these images we identified two point sources, c onsistent 
with the initial XRT position (IVergani et al.ll2010l ). We ac- 
quired spectrosco py of these with the 3-a rm echelle spectro- 
graph X-shooter (|D'Odorico et al.ll2006l ). mounted on Unit 
Telescope 2 of the VLT. These sources were later shown 
to be compact HII regions unrelated to the supernova. In 
this paper we will make use of the spectra obtained of the 
brightest of these tw o regions, referred to as source 'A' in 
IVergani et al.l (|201Cl ) (see Figure [l}. We obtained two spec- 
tra each of 1200 s exposure of source A, using a 5" nod throw 
along the slit to achieve better sky subtraction. We used 
slit widths of 1.0, 1.0 and 0.9" for the UVB, VIS and NIR 
arms, respectively. Data were reduced using version 0.9.4 of 
the ESO X-shooter pipeline ( jGoldoni e t al. 2006). Calibra- 
tion data (bias, dark, arc lamp, flatfield and flexure control 
frames) were taken the same night. Extracted spectra were 
flux calibrated using spectrophotometric standard star ex- 
posures, also taken the same night, resulting in flux call- 
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brated spectra spanning the ultraviolet to the near-infrared 
(~ 300 — 2300 nm) simultaneously. 

We obtained further imaging of the localization with 
Gemini-South on the subsequent three nights, using these we 
were able to identify a fu rther bright point source superim- 
posed on the host galaxy (|Levan et al.ll2O10l ). Image subtrac- 
tion performed between the different images showed a clear 
brightening in this source, which we subseque ntly identi fied 
as the associated rising supernova ((Wiersema et al.ll2010l . see 
Figure [T| . 

We obtained deep imaging of the position of 
GRB 100316D using the FORS2 instrument on Unit Tele- 
scope 1 of the VLT. Seeing conditions during both epochs 
were good, averaging 0.8", but the second epoch observa- 
tions suffer from fairly bright background caused by high 
airmass and a bright moon. The GRB-SN is clearly visible 
in both epochs. 

Following the a nnouncement of a SN discovery 
l|Wiersema et ahlboiOl ). we triggered our cycle 17 HST pro- 
gramme. Here we report the first epoch, obtained on March 
25, ~9 days after the BAT trigger. 

The details of the ground-based and HST observations 
can be found in Table [T] The data were reduced using stan- 
dard procedures in IRAF. 



3 ASTROMETRY AND REDSHIFT 

Subtraction of early from late images taken with Gemini- 
South reveals a variable source inside the X-ray error cir- 
cle which we conclude is GRB 100316D/SN 2010bh (Figure 
[l]). We find an accurate position for this variable source by 
astrometrically calibrating the images, using 190 USNO-B 
stars in the field. We determine the position of the residual 
in the image subtraction by using a Gaussian centroiding. 
This gives a position for the supernova of RA, Dec (J2000) 
= 07h 10m 30.558s, -56d 15' 20.18" (± 0.26"). The error 
is dominated by the scatter of the USNO positions around 
the best fitting polynomial of the astrometric calibration. 
We refine this using an image taken with HST, to our best 
position of RA, Dec (J2000) = 07h 10m 30.530s, -56d 15' 
19.78" (± 0.05"). 

We detect a large number of nebular emission lines 
commonly found in star forming regions with good signal 
to noise in the VLT X-shooter spectrum of source A (see 
Section [T] Figure From the 12 brightest lines we mea- 
sure the redshift z = 0.0591 ± 0.0001. Using the cosmology 
Qm = 0.27, f^A = 0.73 and Ho = 71 km s"^ Mpc"\ the 
luminosity distance is 261 Mpc. 



4 TEMPORAL CHARACTERISTICS OF THE 
PROMPT EMISSION 

The BAT and XRT early emission light curves are shown 
in Figure [2l The 7-ray emission was detected at the IOct 
level at To— 475 s, prior to which there was no detected 
emission at 5a or above. The emission peaked at Tq — 100 s, 
and decayed with an exponential decay constant of 750 s. 
Emission continued through the slew at To-1-750 s. BAT re- 
observed at To -1-5050 s from which time 7-ray emission was 
no longer detected. 
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Figure 2. Top panel: 14-195 keV long-term BAT light curve, 
comprising both event and survey binned data. The zero count- 
rate level is indicated with a dashed line, and the time range 
covered by the panels below is overlaid as a double-sided arrow. 
Lower two panels: zoom in on the 14—195 keV BAT event light 
curve and 0.3-10 keV XRT light curve from To to To+810 s. The 
dashed lines indicate time slices adopted for spectral fitting. 



The prompt emission duration of this GRB is one of 
the longest ever measured. GRB 100316D was detected with 
BAT from ~To-500 s to at least ~To-f800 s, hence the 
lower limit on the duration of GRB 100316D is ~1300 s. 
The fluence, derived from the 15-150 keV spectrum accu- 
mulated between Tq— 475.0 and Tq-I- 795.0 s since trigger, is 
S„ = (5.10 ± 0.39) xlO"*^ erg cm"^ We note that as 7- 
ray emission continued after this time but Swift was unable 
to observe, this should be considered a lower limit to the 
burst fluence. The isotropic equivalent energy that follows 
from this is Ei^o > (3.9±0.3) xlO*" erg. This is valid for 
the energy range 15-150 keV, but if we extrapolate down 
to 1 keV given that this is a soft burst, we find that the 
fluence and therefore Eiso may be higher by 50% leading to 
Eiso > (5.9±0.5) X 10"^ erg. We note that little flux is ob- 
served above 50 keV; the contribution to the isotropic energy 
from very high energy 7-rays is likely to be small. 

The BAT and XRT simultaneous coverage spans 603 s. 
The XRT count rate light curve decays very slowly with 
Q = 0.13 ± 0.03 throughout this interval. 

No variable optical or UV source has been de- 
tected during Swift observations with the UVOT, 
hampered by the bright un de rlying host galaxy 
(|Oates. De Pasauale fc Stamatikoi I2OI0I ). The hght 



curves for each filter are consistent with a constant value, 
and we confirm that no variable source is detected with 
image subtraction using the full UVOT dataset to date. 
Upper limits on the magnitudes of any source at the GRB 
position are given in Table [2l 
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Table 1. Log of Gemini-South, VLT and HST observations. 
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Figure 1. Our Gemini-South discovery images of SN 2010bh, associated with GRB 100316D. The r-band images, obtained one night 
apart, clearly show the afterglow/SN superimposed on its host light. The GRB/SN position, image scale and North and East directions 
are indicated in the middle panel. The star forming region 'A' is indicated in the left panel. The result of PSF matched image subtraction 
is shown in the right hand panel, revealing a brightening of the transient. The Swift XRT error circle is overlaid on this panel. 



5 SPECTRAL CHARACTERISTICS OF THE 
PROMPT EMISSION 

We have carried out time-resolved spectroscopy during the 
593 s of BAT-XRT overlap from To-l-144 s to To-|-737 s. 
For the X-ray data, we take six consecutive time slices, each 
of ~99 s duration. We fit an absorbed power law model, 
and then include a further blackbody component (Figure 
[3]). In either case the X-ray spectral shape appears stable 
throughout the WT observations. The intrinsic absorption 
measured in both sets of fits is consistent at each epoch. 
We then fit the time-averaged WT spectrum consisting of 
data from To -1-144 to To -1-737 s since trigger to obtain the 
best constraints on the intrinsic column which we use in 
the following BAT-XRT joint fits. The time-averaged WT 
spectrum and the time-sliced WT spectra are adequately fit 
with an absorbed power law model, however the inclusion 
of a blackbody provides a significant improvement in the 
fits (F-test probabilities of 10~^° and 10~® respectively), re- 
quiring a larger intrinsic absorption column to accommodate 



Table 2. UVOT 3a- upper limits on the magnitude at the position 
of the SN, from event and image data. We note that these are 
background subtracted but not host-subtracted. The time ranges 
were chosen to overlap with the simultaneous BAT-XRT coverage. 
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Figure 3. Unfolded XRT WT time-averaged spectrum, EF^ vs. 
E, with an absorbed power law plus blackbody model, showing 
the contributions of the two components. 



a kT — 0.14 ± 0.01 keV blackbody component while the 
power law photon index remains approximately the same. 
The blackbody component contributes ~3% of the total 
0.3-10 keV X-ray flux in this model, and has a luminosity 
of (3-4) X 10*^ erg corresponding to an emitting radius 
of 8x10 cm. The blackbody luminosity and temperature 
(Figure [4]) are consistent with remaining constant until at 
least To-HSOO s. 

In order to more rigorously test the significance of the 
blackbody component, we simulated 10000 spectra based 
on the absorbed power law model that best fits the time- 
averaged WT spectrum of GRB 100316D (Tabled). We fit- 
ted the simulated spectra with absorbed power law plus 
blackbody models to recover the chance probability that 
adding a blackbody component significantly improves the 
fit, fo llowing the methodology described in iHurkett et al.l 
l|2008l ). There were no simulated spectra for which Ax^ in- 
dicated an improvement in the fit comparable to that found 
in our data, implying that the blackbody component in the 
time-averaged WT spectrum of GRB 100316D is > 4a sig- 
nificant. This test gives us a high degree of confidence in our 
identification of a thermal component in the X-ray spectrum 
at early times. 

The joint BAT-XRT fits were performed with two 99 s 
and two 198 s consecutive time slices. We fit these spectra 
with an exponentially cut-off power law modeQ with and 
without a blackbody component. A constant normalisation 
ratio between the BAT and XRT spectra (to account for 
calibration uncertainties in the response matrices of the dif- 
ferent instruments) is not included in the reported fits since 
we found this to be 0.9±0.2, i.e. consistent with unity. The 
exponentially cut-off power law model provides a good fit 
with a stable low energy power law slope of F ~ 1.4. The 
peak energy is at -^30 keV in the first time interval and 
moves to lower energies with time (Figure |4)). To estimate 
the spectral shape during and prior to the BAT trigger we 
fitted the BAT spectrum from To-175 to To-|-144 s simuF 
taneously with the time-averaged WT spectrum, given that 
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Figure 4. Evolution of power law photon index F, peak en- 
ergy iJpit and blackbody temperature kT in the blackbody plus 
exponentially cut-off power law model fitted to the BAT-XRT 
spectra for GRB 100316D (solid error bars), compared with the 
same model fitted to BAT-XRT data of GRB 060218 (taken from 
Kaneko et al. 2007, dotted error bars) in the source rest-frames. 
Errors are 90% confidence. The bottom panel shows the percent- 
age contribution of the blackbody component to the total ob- 
served 0.3—10 keV X-ray flux, in time-sliced XRT spectral fits to 
GRB 100316D. 



1 defined as A{E) = 



no significant X-ray spectral evolution is observed in the 
interval 144-737 s. We stress that because the data are not 
simultaneous this fit provides only an indication of the spec- 
tral shape at early times. Spectral fitting results are given 
in Table O 

The latest time slice taken for our joint BAT-XRT fits 
also has full coverage in the u band with UVOT. We have 
included the u band 3a upper limit and refit the data with 
an absorbed and reddened blackbody-|-power law model. We 
obtain an acceptable fit, when the underlying spectrum is 
fixed in shape as listed in Table |3l with an intrinsic extinc- 
tion of E{B — V) ^ 0.9 mag or Av ^ 2.6 mag using the 
Small Magellanic Cloud extinction law (|Peil il992. fixing the 
Galactic extinction to E(B — V) = 0.12 mag according to 
ISchleeel. Finkbeiner fc DavisI [l99l l . This provides a lower 
limit on the host galaxy extinction along the line-of-sight to 
the GRB, only under the assumption that any optical tran- 
sient emission lies on an extrapolation of the high energy 
spectrum. A measurement of the optical extinction along 
the line-of-sight to the GRB can be made with spectroscopy 
of the supernova, for example through the strength of inter- 
stellar medium Na I lines (Bufano et al. in preparation). 
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Table 3. Spectral fits to the Swift BAT and XRT data with an absorbed power law (PL) or cut-ofi power law (CPL) model, with and 
without a blackbody (BB) component. Errors are quoted at the 90% confidence level. Galactic absorption is fixed at TxlO^" cm~2. 

* These parameters are tied together for all time intervals. 

* XRT data are not available during this time. We perform this fit only to provide an indication of the cut-off power law parameters over 
this time interval, using the time-averaged WT spectrum given that no significant X-ray spectral evolution is detected in the interval 
144-736 s. Therefore this fit was performed separately rather than simultaneously with the 144-736 s interval BAT-XRT time-sliced fits. 





PL 


or CPL 




PL-f BB or CPL-I-BB 




T 

(s since BAT trigger) 


A^H.int 

(10^2 cm'2) 


r 


(keV) 


A^H.int kT T 
(10^2 cm-2) (keV) 


(keV) 


XRT WT time-averaged: 
144-737 


0.91±0.05 


= 587.9/480 
1 40+"-''4 




X^/v = 538.8/478 
1.83±0.06 0.14±0.01 1.53±0.06 





XRT WT time-sliced: 






= 830.9/754 




X^/y = 782.1/742 


144-243 


0.90±0.05' 


1.35±0.06 - 


1.8±0.2« 


0.14±0.02 


1.5±0.1 


243-342 


II 




1.41±0.06 - 


// 


0.14±0.02 


1.6±0.1 


342-440 


n 




1.36±0.06 - 


// 


0.15±0.02 


1.5±0.1 


440-539 


II 




1.43±0.06 - 


// 


14+0'03 
'^■^^-0.02 


1.6±0.1 


539-638 


II 




1.43±0.06 - 


// 


0.13±0.02 


1.5±0.1 


638-737 


II 




1.40±0.06 - 


// 


14+003 
'^■^^-0.02 


1.6±0.1 



XRT PC time-averaged: 
32906-508263 


C-statistic/i/ = 13.8/15 
0.90 fixed 3.5±1.0 


C-statistic/;^ = 15.4/15 
1.83 fixed - 4.6l^-^ 


BAT: 
-175-144 


X^/y = 16.0/14 
2.1±0.4 





BAT-XRT: 


x^/y 


= 617.4/518 


3311 




X^/y = 572.1/516 




-175-144* 


0.90 fixed 


1.33±0.03 


1.83 fixed 


0.14±0.01 


1 45+0.02 


BAT-XRT: 


X^/y 


= 853.4/768 


25+8 




X^/y = 806.1/760 


28j;f 


144-243 


0.90 fixed 


1.25±0.08 


1.83 fixed 


0.14±0.02 


1.4±0.1 


243-342 


// 


1.32±0.08 


22t^ 


// 


0.14±0.02 


1.5±0.1 


26+f 


342-539 


// 


1 24+0-06 


14ti 


// 


0.15±0.01 


1.3±0.1 


13+^ 


539-737 


// 


1.28±0.07 


14±3 


// 


0.14±0.01 


1.4±0.1 


13^^ 



6 THE LATE-TIME BEHAVIOUR 

The X-ray emission after 10'' s decayed more steeply than 
did the emission during the first 800 s after the BAT trigger. 
Treating the early and late X-ray emission separately, we fit 
only the data after 10'* s with a power law decay. This results 
in an acceptable fit, with a = 1.35 ± 0.30. 

The XRT count rates for this source lie within the 
spread of the entire XRT GRB light curve sample at all ob- 
served times, however its shape is atypical. Four light curve 
types were identified in the Swift XRT catalogue, to whic h 
76% of their sample could be ascribed (lEvans et al.ll2009h . 
GRB 100316D appears not to fit any of these light curve 
types, if all the early data are included. 

The X-ray spectrum from To-|-3xlO'' to SxlO'' s al- 
lows only poor constraints on the parameters of an ab- 
sorbed power law model, but we find it is significantly spec- 
trally softer at this time than in the earlier WT data. Us- 
ing Cash statistics, and fixing the intrinsic absorption to 
that found in the combined "WT spectrum, we find the 
power law has a photon index of F = 3.5 ± 1.0 (Table |3]). 
This is a very high value, or a very soft spectr um, that al- 
thoug h not unique among Swift bursts (see e.g. lEvans et al.l 
I2OO9I ) is difficult to accommodate within the standard fire- 
ball model generally used to expla in GRB afterglow emission 
(e.g. ISari. Piran fc Naravan|[l998l ), but may possibly be ex- 
plained with inverse compton scattering (e.g. Waxman et al. 



2007). The temporal decay rate, however, is rather typical 
of GRB afterglows. 

We now turn to the optical regime, to estimate the su- 
pernova contribution at the epochs of our VLT data. We 
performed point spread function (PSF) matched image sub- 
traction on the F0RS2 images, using the ISIS2 code (Alard 
2000). For the second epoch Vhigh data we use the short 
exposure time images, as the frames with longer exposure 
time show the supernova to be somewhat saturated. The 
subtraction shows a residual at the position of the supernova 
associated with this GRB l|Wiersema et al.ll2010l ) in the V 
and I bands, whereas the subtraction is too noisy in the R 
band. This residual shows that the supernova is fainter in the 
second F0RS2 epoch than in the first epoch, at least in V 
and I, suggesting that the supernova peaked in these bands 
between 8 and 13 days after burst. This is similar to the 
SN2006aj (GRB 60218), which peak e d at 10.2 ± 0.3 days 
after burst in V (iFerrero et all feooel: iMirabal et all I2OO6I : 



Modiaz et aLll2006l : ISollerman et al.ll2006l ) 



T THE HOST GALAXY 

The deep Gemini, VLT and HST observations show a bright, 
blue galaxy, covering the 3.7" radius refined XRT error circle 
(Figures [T] and |5| . The galaxy has a highly disturbed mor- 
phology, but is dominated by a central bright region, pos- 
sibly the galaxy nucleus, and an arc-like shape resembling 
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a spiral arm. The galaxy hosts a large number of knots, or 
blobs, likely giant HII regions. Colour differences between 
these knots can be discerned, caused by differences in age, 
metallicity, reddening or a combination of these. The super- 
nova itself is located on top of a bright, blue knot, close 
(in projection) to the possible galaxy nucleus. The UVOT 
also observed the host galaxy in all filters and approximate 
magnitudes, uncorrected for extinction and including an un- 
known level of contamination by the afterglow and/or SN, 
are reported in Table ID 

It has long been contended that the properties of the 
host galaxies of GRBs can provide additional clues to the 
progenitor evolution, through determination of properties 
such as the age, metallicity and alpha element enrichment 
of the stellar population. The low redshift GRB-SNe are 
of particular interest, as they allow one to obtain spatially 
resolved spectroscopy (e.g. ^Christensen et al. 2008). The de- 
rived properties can then be meaningfully co mpared to those 
of 'normal' core-collapse supernova hosts (|Fruchter et al.l 
l2006l : [Svensson et al.ll2O10l ) or other star forming galaxy pop- 
ulations. For these reasons we analyse the X-shooter spec- 
trum we obtained of a bright HII region w ithin the host, 
which we call source A (|Vergani et al.ll2010l ). spatially dis- 
tinct from the supernova region (Figures[T]and[8l). We choose 
to use this spectrum rather than the X-shooter spectra of 
the supernova itself, as in these spectra the contribution of 
the supernova to the spectrum is large and can not be easily 
separated from the underlying HII region lines. Source A 
is located near the supernova (projected distance approxi- 
mately 2kpc) and we may expect the properties of this re- 
gion to give a meaningful insight into the properties of the 
host as a whole. 

We detect a large number of nebular emission lines com- 
monly found in star forming regions; a selection of those 
used here (which are all unresolved) is shown in Figure [S] 
We use the Balmer Ha and H/3 line fluxes to measure a value 
for the Balmer decrement (jointly fitting the stellar Balmer 
absorption component and the nebular emissio n line com- 
ponent), assuming case B r ecombination (e.g. lOsterbrockl 
ll989l : llzotov. Thuan fc Lipov etskv 1994). Assuming an elec- 
tron temperature of 10'* K we compute the extinction (we 
find E{B — V) = 0.178 mag for the combination of Galactic 
reddening of E{B — V)gax = 0.07 mag and intrinsic reddening 
in source A), and deredden the spectrum accordingly. 

We compute the average electron density Ue using the 
[SII] doublet and the [Oil] doublet: both of the doublet 
fiux ratios are consistent with their lower li mits, so we 
assum e a value of = 100 cm~^ (see e.g. lOsterbrockl 
Il989l ). The electron temperature sensitive [O III] A4363 and 
[O II] AA7320,7331 doublet are detected, from which we cal- 
cula te the electron te mperatures for O^"^ and O"*" directly 
fe.g. lOsterbr"oc^ll989l ). at 11880±800 K and 10400±1100 K, 
respectively (errors are lir). lonisation corrections to ac- 
count for the presence of O'^^ were found to be negligi- 
ble, also evidenced by the non-detection of Hell A4686, 
so we determine the oxygen abundance in region A as 
0/H = 0+/H+ + 02+/H+. Using the measured line fluxes, 
densities and temperatures we arrive at an electron tempera- 
ture based oxygen abundance of 12 + log(0/H) = 8.23±0.15 
in region A (0.37 times Solar metallicity, assuming a so- 
lar oxygen abundance of 12 4- log(0/H) = 8.66; Asplund 
et al. 2004). We note that a similar oxygen abundance has 



Table 4. UVOT observed magnitudes and la errors for the host 
galaxy. Photometry is performed as described in Section [2.1l N.B. 
any afterglow or supernova contribution is uncertain so could not 
be removed, however no variability is seen in subtraction of late- 
time images from earlier epochs. 
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been derived from spectra at the SN site fsee lChornock et al.l 
I2OIQ) . 

In the spectrum we detect Balmer absorption lines un- 
derlying the nebular emission. This is common in GRB host 
galaxy spectra, and can provide a useful age tracer, pro- 
viding the spectra have sufficient resolution (e.g. Wiersema 
et al. in preparation). We fit the H(5 absorption compo- 
nent, and find an approximate age for the dominant stel- 
lar popula tion of ~ 30 Myr, assuming contirmous star 
formation (|Gonzalez Delgado. Leitherer fc H cckman 1999). 
We also note the absence of obvious Wolf-Rayet star emis- 
sion line features in source A. The detection of the bright Ha 
line allows us to estimate the star formation rate at location 
A: the Ha line luminosity is a tracer of recent star formation 
frequently used in nearby GRB host galaxies. We use the ex- 
pression SFRna = 7.9 x 10"*^Lhc« (|Kennicuttll 19981 ). Using 
the de-reddened flux from the X-shooter spectrum we flnd 
SFRua = 0.17 Mq yr~^. Note that no aperture correction 
is performed, so this value is formally a lower limit. 

A more detailed analysis of the host galaxy will appear 
in Flores et al. (in preparation). 



8 DISCUSSION 

GRB 100316D is an atypical gamma-ray burst both in its 
temporal and spectral behaviour. The very soft spectral 
peak and extended and slowly decaying flux emission are 
highly un usual amon g the prompt emission of GRBs (e.g. 
iSakamoto et al.ll2008l ). The estimated total isotropic energy 
(£iao > (5.9±0.5) xlO*^ erg), when considered with the 90% 
confldence range of spectral peak values of 10-42 keV we 
measure in Section [5] places GRB 100316D on the low i?iso- 
low i5pk tail of the Amati correlation wh ich relates -Ep k and 
Eiso for the majority of typical GRBs (|Amatil I2OO6I ). The 
nota ble outher to this relation is GRB 9804 25 (SN 1998bw) 
(e.g. iKaneko et"ai1l2007l : I Amati et aLll2007h . We also com- 
pare t his GRB with the relation reported bv lGehrels et al.] 
(|2008l ) between XRT flux at 11 h (here, (3.2±0.8) x 10"*^ erg 
cm~'^ s~*) and BAT fluence and flnd, using the lower limit 
to the fluence, that it lies along the correlation with values 
that are typ i cal fo r the sample of Swift GRBs reported in 
lEvans et all (|2009l ). 

GRB 100316D is among the longest duration Swift 
GRBs (Tgo 1300 s), after GRB090417B (T90 ~ 2130 s) 
and GRB 060218 (T90 - 2100 s). GRBs with durations of 
400 s or longer, as quantifled using the T90 parameter, ac- 
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Figure 5. A selection of emission lines wliich are most im- 
portant for the oxygen abundance analysis, detected in tlie X- 
shooter spectrum of the source 'A' HII region in the host galaxy, 
~2kpc projection from the GRB/SN location. The top two pan- 
els show the Balmer Ha and H/3 lines used for the determination 
of the reddening. The middle panels show the strong lines of the 
[Oil] AA3727,3729 doublet and [O III] A5008. The zero velocity 
point for [O II] AA3727,3729 is (arbitrarily) chosen to be the cen- 
tre of the redmost component of the doublet. The lower panel 
shows the temperature sensitive lines of [O III] and [O II] (of the 
latter only one member of the doublet is shown) from which we 
derive T^{0'^+) and Te(0+) respectively. 



count for only 1.4% of the Swift sample observed to date. 
Owing to their long prompt emission durations, these GRBs 
are more likely to have exceptional broadband coverage, 
lending themselves to prompt emission studies. We show 
the X-ray light curve and hardness ratio of GRB 100316D 
in relation to the Tgo > 400 s XRT-observed Swift events 
in Figure [5] (omitting GRB 060218 as it can be seen in Fig- 
ure [T] data are taken from the Swift XRT GRB light curve 
repositor>0, lEvans et al]|2007l ). 
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Figure 6. Upper panel: X-ray light curves in flux density at 
1.7 keV against time since trigger for a selection of Swift XRT- 
observed very long prompt emission duration (>400 s) GRBs (ex- 
cluding 060218). Lower panel: (1.5-10 keV)/(0.3-1.5 keV) hard- 
ness ratios derived from the XRT count rate light curves. 



This category comprises a number of different light 
curve types, from those which are long due to triggering on 
a precursor or the presence of a broad flare to the very grad- 
ually decaying prompt emission of GRB 100316D. The Swift 
subsample shown here indicates two GRBs, in addition to 
GRB 060218, whose emission resembles GRB 100316D both 
in temporal decay and in spectral hardness: G RBs 070616 
( Starling et al. 2008) and 090419 (|Stratta et al., 2009). We 
note that GRB 090417B is also spectrally hard at early 
times, but this may be attributed to strong X-ray flaring 
(Holland et al. in preparation). 

GRB 070616 was particularly well sampled, with 7-ray 
through optical coverage at early times clearly showing the 
movement of the spectral peak energy of the Band function 
(|Band et a l."l993') and the occurrence of additional spectral 
softening (Starling ct al. 2008). However, the spectrum was 
much harder than that of 100316D - observed by both Swift 
BAT and Suzaku WAM, and neither the redshift nor the 
presence of any SN could be established as the source lay 
close to a bright star. 

We now discuss the properties of GRB 100316D in 
the context of GRB-SNe, paying particular attention to 
GRB 060218. 



8.1 Comparison to GRBs with accompanying 
supernovae: high energy emission 

In Figure [7] we show the X-ray light curves and hardness 
ratios of a subset of the GRBs with associated SNe (ei- 
ther spectroscopically or photometrically confirmed). On 
this plot we also include GRB 051109B which w as tentatively 
assoc iated with a very low redshift (2 = 0.08, IPerlev et al.l 
l20()f|) host galaxy and has no reported SN that we are 
aware of, and the well-studied nearby GRBs 060505 and 
060614 for whic h SN searches to very deep limits registered 
no detections (|Fvnbo et al.l l2006l : lOella Valle et al.l I2OO6I : 
iGal-Yam et al.ll2006l ). The prompt emission shows a broader 
range of observed fluxes than the late-time decays, perhaps 
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suggesting they are unrelated or that an extra component 
is present in the early emission which contributes differently 
in each GRB. The majority of the GRB-SNe (and GRBs 
without SNe shown here) light curves decay steeply over the 
first 1000 s, with GRBs 100316D and 060218 being the only 
exceptions. The hardness ratios reveal that GRB 100316D 
and GRB 060218 have a similar spectral hardness evolution, 
clearly different from those of most other GRB-SNe which 
remain at a stable and softer spectral shape throughout the 
prompt and late-time emissions. GRBs 060614 (without a 
SN) and 090618 (photometrically discovered SN) also tran- 
sition from hard to soft, but do this earlier, over the first 
200 s during their steep decays. We note, however, that 
X-ray data during the early emission for a number of the 
GRB-SNe, including the spectroscopically confirmed GRB- 
SNe 980425, 030329, 031203 and 050525A, are either not 
available or inadequate for spectral studies. 

8.2 Comparison to GRBs with accompanying 
supernovae: host galaxies 

Nearby GRBs show a large variety in their host galaxy prop- 
erties: the prototype GRB-SN, GRB 980425 and SN 1998bw, 
occured in a dwarf spiral, in a small H II region right next to 
a very large star formation complex; whereas GRB 060218 
and SN 2006a,i, occured in a v ery faint, blue compact dwarf 
galax;y l|Wiersema et al.ll2007^ . This diversity is illustrated 
in Figure [8] 

Assuming that the properties of source A are at least 
roughly representative of the host galaxy as a whole, we can 
compare the spectroscopic properties with those of the other 
nearby GRB-SN hosts. Electron temperature based oxygen 
abundances have been derived for several of these, see Table 
[5] The case of GRB 980425 is of particular interest: besides 
many similarities between SN 1998bw and SN 2010bh there 
are also several similarities in host galaxy properties. The 
metallicity measured at the location of the GRB and at the 
locat ion of a nearby, bright, WR star region l|Hammer et al.l 
I2OO6I ) are clearly similar to the metallicity measured for the 
HIT region A in the host of GRB 100316D. In addition, the 
size and brightness (and possibly morphology) of these two 
hosts are comparable, though an important difference be- 
tween the host of 980425 and 100316D is evident from Figure 
[8l the morphology of the host of 100316D is highly disturbed, 
possibly indicative of a recent merger. 

The host of GR B 980425 has been studied using integral 
field spectroscopy (|Christensen et al] I2OO8I ) , showing that 
the HII region in which the GRB exploded is similar in 
mass, SFR, reddening and line equivalent width to other 
H II regions in this host, with the exception of a bright WR- 
star rich H II region 800 pc away from the GRB site. In 
terms of metallicity, the WR region and GRB site have a 
some what lower metallicity t han other HII regions in this 
host (jChristensen et al]|2008h . If this situation is compara- 
ble to the host of GRB 100316D, we may expect that the 
properties measured above of source A are indeed represen- 
tative of most of the host galaxy properties. 

8.3 Comparison to GRB 060218 

In Section 18.11 we found that among the GRB-SNe, there 
is a clear commonality in the X-ray behaviour of 100316D 



and 060218. Comparing the prompt emission spectra in Sec- 
tion [5] we again find parallels in that both bursts seem to 
require a similar thermal component in addition to the typ- 
ical GRB synchrotron emission for which the synchrotron 
peak is observed to move to lower energies with time. The 
time-evolving prompt spectral parameters are compared in 
Figure |4] Both events have low isotropic equivalent energies, 
of order 4x10'*^ erg, and very little flux above 50 keV. 

These two events lie in apparently qui te different envi- 
ronments: 060218 in a fain t, compact dwarf (jWiersema et al.l 
I2OO7I : iModiaz et ahl 12006*) and 100316D in a luminous, dis- 
turbed possibly spiral host galaxy. The metallicities are also 
different, with the host of 100316D likely being more metal- 
rich (though the metallicity is below Solar) . A prompt opti- 
cal component, as observed in 060218, could not be detected 
in 100316D, but conditions for its detection were far less op- 
timal in this case due to the superposition on a brighter 
host galaxy and the possibility of a higher extinction along 
the line-of-sight. The higher intrinsic X-ray column density 
we measure for 100316D compared with similar modelling of 
0602 18 does not necessarily imply a higher optical extinction 
( Watson et al.ll2007l : [Campana et al.ll2010l ). 

In summary, both 100316D and its predecessor 
060218 are nearby {z = 0.059,0.033), long-duration 
(Tgo ^1300 s, 2100 s), initially relatively constant in X- 
ray flux (Figure |6]), spectrally soft (very few or no counts 
above 100 keV; low -Epk, Figure subenergetic (both have 
Eiso ~ 4 X lO''^ erg) GRBs with an associated SN. These 
two events show similar prompt emission properties and 
stand out among the GRB-SNe subsample considered here 
for their unusual X-ray evolution. However, their host galax- 
ies are rather different in morphology and metallicity, with 
the host of 100316D more closely resembling the host of 
980425. 

The thermal X-ray component, with a luminosity, tem- 
perature and radius similar to that observed in 060218, dic- 
tates that the shock break-out of the supernova must be 
considered. The optical/UV thermal component observed 
in 060218 is not seen here: if the extinction is similar we 
would expect a shock breakout similar to 060218 to be two 
or more times fainter in 100316D and with the relatively 
brighter host galaxy it is perhaps no surprise that early op- 
tical emission is not detected with UVOT. The presence of a 
thermal component in the prompt 7-ray emission of BATSE 
(Burst And Transient Source Experiment) G RBs, in add i- 
tion to non-thermal emission, was proposed bv lRvdd l|2004h . 
The discovery of thermal components in th e Swift XRT X- 
ray spectra shown here for 100316D and in I Campana et al] 
(2006) for 060218 will therefore be important to consider in 
future studies of GRB prompt emission. The dominant com- 
ponent of the high energy emission in 100316D remains the 
synchrotron-like non-thermal spectrum common to all types 
of GRB (with and without SNe) thought to originate in in- 
ternal shocks in a relativistic jet. The long duration of the 
early X-ray emission is curious, and exceedingly rare, per- 
haps suggesting a greater reservoir of material is available 
to feed the central engine and prolong its activity. 

The discovery of GRB 100316D and its associated su- 
pernova SN 2010bh and the analyses of its high energy emis- 
sion and host galaxy properties presented here illustrates the 
diversity in GRB-SNe characteristics that must be under- 
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Figure 7. X-ray light curves against time since trigger in fiux density at 1.7 keV for a selection of long GRBs with associated supernovae 
(either spectroscopically confirmed - generally coloured in shades of red with GRB 100316D in black, or photometrically determined - 
shades of blue and green), two nearby long GRBs with no observed supernovae to deep limits and a further long GRB likely at z 0.1 
with no reported supernova (grey shades). The lower panel shows the hardness ratios for the Swift XRT-observed GRBs, derived from 

the 1.5-10 keV and 0.3-1.5 keV count rate light curve s in the same col our scheme^ 

Reference s for the data sour ces: lOOSieD^, 980425 llPian et al.ll2000l). 030329 llWillingale et al.ll2004). 03 1203 l lWatson et al. 
050525A llBlustin et al.ll2006l. and 2), 060218^, 081007^, 970228 jCosta et al.lll997l : jPe Pasguale et al.ll2006l) . 011121 ijPiro et al. 
0508242, 070419A2, 080319B2, 090618^ 051109B2, 060505^, 0606142. 



2004), 



2005) 



Table 5. Host galaxies of nearby GRB-SNe. The absolute magnitudes Mg were taken from the compilation of Levesque et al. (2010). 



GRB Te oxygen abundance Host type Absolute magnitude References 

(12 + log(0/H)) Mb 



980425 


8.25 (GRB site) 

8.39 (nearby WR region) 


Dwarf spiral 


-17.6 


SoUerman et al. (2005); 
Hammer et al. (2006) 


020903 


7.97 


Irr 


-18.8 


Hammer et al. (2006) 


030329 


7.72 


Irr 


-16.5 


Levesque et al. (2010) 


031203 


8.02 ± 0.15 (integrated) 


Irr 


-21.0 


Prochaska et al. (2004) 


060218 


7.54l°'j^ (integrated) 


Irr 


-15.9 


Wiersema et al. (2007) 


100316D 


8.23 ±0.15 (source A) 


Spiral? Irr? 


-18.8 


This work 



stood if we are to fully appreciate the relationship between 
GRBs and core-collapse SNe. 
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